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ABSTRACT 

Aims. We present the results of the X-ray spectral analysis of an XMM-Newton-selected type II QSO sample with z > 0.5 and 0.5- 
lOkeV flux of 0.3 - 33 x 10~ 14 erg/s/cm 2 . The distribution of absorbing column densities in type II QSOs is investigated and the 
dependence of absorption on X-ray luminosity and redshift is studied. 
Methods. We inspected 51 spectroscopically classified type II QSO candidates from the XMM-Newton Marano field survey, the 
XMM-Newton-2dF wide angle survey (XWAS), and the AXIS survey to set-up a well-defined sample with secure optical type II iden- 
tifications. Fourteen type II QSOs were classified and an X-ray spectral analysis performed. Since most of our sources have only ~40 
X-ray counts (PN-detector), we carefully studied the fit results of the simulated X-ray spectra as a function of fit statistic and binning 
method. We determined that fitting the spectra with the Cash-statistic and a binning of minimum one count per bin recovers the input 
values of the simulated X-ray spectra best. Above 100 PN counts, the free fits of the spectrum's slope and absorbing hydrogen column 
density are reliable. 

Results. We find only moderate absorption (Nn = (2 - 10) x 10 22 cirT 2 ) and no obvious trends with redshift and intrinsic X-ray lumi- 
nosity. In a few cases a Compton-thick absorber cannot be excluded. Two type II objects with no X-ray absorption were discovered. 
We find no evidence for an intrinsic separation between type II AGN and high X-ray luminosity type II QSO in terms of absorption. 
The stacked X-ray spectrum of our 14 type II QSOs shows no iron Key line. In contrast, the stack of the 8 type II AGN reveals a very 
prominent iron Ka line at an energy of ~ 6.6 keV and an EW ~ 2 keV. 
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<0 1 . Introduction thick torus. Furthermore, their lack of emission lines in a wide 

25 ■ optical wavelength range hampers their identification. About 

O ■ The unification model for active galaxy nuclei (AGN) is based 150 optically selected type n QSOs from the sloan Digital 

L; ; on the assumption that the Seyfert type I/type II dichotomy in sky Suryey wer£ studied in Zakamska ^i] J200I . However 

. _H . AGN is a result of varying orientation relative to the line of sight type u QSO s are more efficientl y found in follow-up observa- 

^ ■ of similar objects (| Antonucci|1993|) . All AGN consist of a super- tionso f X -ray surveysJSzokolv et al. 2004; Mai nieri et alJ 2002: 

H ; massive black hole with an accretion disk. This central engine is iKrumpe et al.ll2lMlBarcons et al]l20071. iTedds etdF oOS). 

^ surrounded by optically thick, toroidally concentrated dusty ma- 

terial. An observer can either view the central engine directly T yP e 11 Q SOs P la ? an important role in understanding the 

(type I AGN) or through the optically thick torus (type II AGN). X " ra y universe. Since they show significant absorption, type II 

Direct observation of the broad line region reveals a strong UV QSOs/AGN are considered to be a main contributor to the hard 

continuum and broad permitted emission lines while the optical X " ra y background and their existence in considerable numbers 

spectra of type II AGN show narrow permitted and forbidden iLDeeded for the synthesis of the X-ray background d Gilli et al. | 
emission lines. lHao et a l. (2005) provided a formal separation Posslble evolution of the absorption with intrinsic lurm- 

criterion by finding a bimodal disttibution in the Ha-FWHM noslt y and / or redshlft 18 a matter of intensive debate and would 

wit h a significant dip at FWHM(Ha)=1200km/s. essentially influence the X-ray background synthesis models. 

iKleinmann et all (1 19881) discovered the first type II quasi Many papers discuss the fraction of absorbed AGN as a 

stellar object (QSO). QSOs show higher intrinsic luminosity function of luminosit y and/or redshift, and this is a highly 

than AGN. In the optical the conventional dividing line between controversial subje ct dDwellv & Pagd 120061 : iTozzi et alj [2006; 

Seyfert galaxies and QSOs is M B = -23. Despite their high in- Akvl as et ai1l2006l). H owever, there are only a few studies (e.g. 

ttinsic luminosities, type II QSOs are very hard to identify. A sig- lLa Franca e t alj 120051) which investigate the evolution of ab- 

nificant fraction of the emitted power is absorbed by the optically sorbing column densities in type II AGN up to high redshift. 

Other studies focus on the local universe, where the existence 



Send offprint requests to: M. Krumpe of so-called Compton-thick absorbed ob j ects is well know n 

* Based on observations obtained with XMM-Newton, an ESA sci- dMaiolino etalj|1998|; |Bassani et al]|1999|; |Risaliti et al.||2000|). 

ence mission with instruments and contributions directly funded by A Compton-thick object is absorbed by column densities of 

ESA Member States and NASA. Nu > 1.5 X 10 24 crrT 2 , so that the cross section of Compton 
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scattering overcomes the photo-electric absorption. Hence, the 
reflected X-ray component is observable. 

There is evidence for a large fraction of Compton-thick ab - 
sorbed quasars at high redshifts (Martinez-Sansi gre et al.l l2007). 
but only one candida te has been reported so far at high redshift 
dNorman et al.l [2002b . and even in this case the column den- 
sity is quite uncertain. Most of the X-ray sources from deep 
Chandra and XMM-Newton surveys are absorbed by column 
densities 10 21 < Na/cm~ 2 < 10 23 , with hard X-ray luminosities 
of 10 42 < L x „„ c /(erg/s) < 10 44 (IComastrill2004t ISzokolv et al l 
2004; lMainieri et al.ll2002l;lMateos et al.ll2005l) . These values are 
also found in Compton-thin absorbed (10 21 < Nw/cmT 2 < 10 24 ) 
Seyfert II galaxies at low redshifts. 

The definiti o n of type II QSOs is somewhat arbitrary. 
Zakamsk a et al.l d2003b define a type II QSO only based on 
its optical properties. They select AGN with narrow permitted 
emission lines (FWHM < 2000 km/s) and classify objects with 
[Om] /15008 line lum inosity of L > 3 x 10 8 L as type II QSO. 
Main ieri et al.l (120021) introduce the term "type II QSO region". 
Although their type I/II classification is based on optical spec- 
tra, a type II QSO has to have an intrinsic X-ray luminosity 
of Lx mT > 10 44 erg/s (0.5-10keV) and an absorbing hydrogen 
column density of Nh > 10 22 crrT 2 . Throughout the paper we 
use the definition of a type II QSO to be an object with nar- 
row forbidden and narrow permitted emission lines in the opti- 
cal spectrum, as well as a de-absorbed intrinsic X-ray luminosity 
L XlNT > lO^erg/s (0.5-10keV band). However, the Seyfert type 
II classification of our sample objects is strictly based on the op- 
tical spectra. 

The present work focuses on properties of the high X-ray lu- 
minosity type II QSOs found in medium deep surveys performed 
by us. The paper is organised as follows. In Sect. [2] we describe 
how we selected our type II QSO candidate sample and we study 
the properties in Sect. [3] The X-ray data and the extraction of 
the X-ray spectra are described in Sect. |4] Since X-ray spectra 
with small numbers of counts have to be fitted, we studied dif- 
ferent fitting and binning methods in Sect. [5] and then describe 
the analysis of the X-ray spectra. The results are discussed in 
Sect. [6] Finally, our conclusions are outlined in Sect. [7] 

Unless mentioned otherwise, all errors refer to a 68% con- 
fidence interval. We assume Hq = 70km/s/Mpc, Qm=0.3, and 
Q A =0.7. 



2. Definition of the type II QSO sample 

We began with a list of 51 spectroscopically classified type II 
QSO candidates which are associated with X-ray sources. All 
optical counterparts have z > 0.5 and observed X-ray lumi- 
nosities (not corrected for intrinsic absorption) of log(Lx OB s) > 
42.5. The type II QSO candidates were taken from three dif- 
ferent X-ray survey s, the XMM-Newton Marano field survey 
dKrumpe et alj|2007l) . the XMM-2d¥ wide angle survey (XWAS, 
Tedd s et al.l 120081). and from 'An XMM-Newton Internationa l 
Survey' (AXIS. iBarcons et al.ll200l 120071: ICarrera et al.ll2007h . 
To obtain a well-defined sample with secure optical type II 
identifications and X-ray counterparts we applied a two stage 
process. First, we visually inspected the optical spectra of all 
type II QSO candidates from all X-ray surveys and determined 
the FWHM of the emission lines. The type II QSO candi- 
date classification was only based on the optical spectra. The 
spectral resolution of the optical setup was ~ 1050 km/s for the 
Marano field survey, -700 km/s for XWAS, and -300-600 km/s 
for AXIS. Therefore, an intrinsic line of 1200 km/s was observed 



with -1250-1600 km/s. We defined emission lines with a mea- 
sured resolution (including instrumental resolution) of FWHM 
< 1500 km/s as narrow. 

The optical spectra were only used to determine a secure 
Seyfert type of the objects. They were not used to separate be- 
tween AGN and QSOs. None of the 51 type II QSO candi- 
dates showed obvious broad emission lines in the optical spectra. 
However, since we aimed for a very strict type II QSO sample, 
we excluded many objects that showed only the low excitation 
On line. The optical spectrum had to comprise several emission 
lines that allowed us to establish an optically secure classifica- 
tion of the Seyfert type. In addition the signal-to-noise ratio had 
to be appropriate in order to verify the non-existence of broad 
permitted emission line. We excluded all doubtful cases. 

We introduced two different categories of type II QSO candi- 
dates based on the reliability of the type II classification. Objects 
that have optical spectra with high excitation lines and at least 
one permitted AGN emission line that is detected but narrow 
(e.g. Ly-a, Civ, H/?) are marked with the optical flag 'S' (see 
Table[TJ - signifying a secure optical identification. Their optical 
classification is very robust. 

Objects with less secure type II identification belong to the 
tentative sample (optical flag 'T'). The classification as tentative 
object can be due to the following reasons: 

- The optical spectrum shows narrow emission lines, com- 
mon AGN high excitation lines, but no Mg n emission. Only 
[Ne v] and [O n] are detected in the spectra. The presence of 
a common AGN high excitation line ([Ne v]) without a broad 
Mg ii suggests the classification as a type II QSO candidate. 
However, since no other permitted emission line is accessi- 
ble, a secure type II classification cannot be established. 

- Common AGN high excitation lines are narrow but show in- 
dications of underlying broad components (optical spectra in 
AppendixlAl. These objects are likely to be transition objects 
between Seyfert type II / type I 

- The signal-to-noise ratio of the H/3 line, or its coincidence 
with the atmospheric band prevents a robust determination of 
the FWHM and could hide a weak, but broad H/? component. 

For all sources in the tentative sample we give an individual 
comment at the end of this section. After inspection of the op- 
tical spectra we were left with 14 secure objects and 12 type II 
objects that belong to the tentative sample. 

As a second step to setting up a well-defined type II QSO 
sample, we verified that the spectroscopically identified counter- 
part is associated with the X-ray source. For a significant num- 
ber of objects much deeper imaging data have become available 
subsequent to the epoch when the spectra were obtained. For ex- 
ample, the XWAS optical counterparts were origi nally selected 
from SuperCosmos optical imaging survey data dHamblv et al.l 
2001 ). For the majority of the objects we now have additional, 
deeper imaging data in different bands. These imaging data were 
obtained with the Wide Field Imager (WFI, 2.2 m telescope at La 
Silla) or the Wide Field Camera (WFC, 2.5 m Isaac Newton tele- 
scope on La Palma). We visually investigated the best available 
imaging data (WFI, WFC, SuperCosmos) for additional optical 
counterparts and rejected doubtful counterpart identifications. 
Furthermore, we computed the probability that the optical ob- 
ject with the /?-band magnitude tor is associated with the X-ray 
source. This was based on 

O exp(-distn V /2) 
2no-tN(< tor) 
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where distox is the distance between optical counterpart and 
X-ray source and <x x is the p osition error of the X-ray source 
(Sutherland & Saun ders! 1 1992b . A detailed description is given 
in iRrumpe et alj d2007l) . Most of our optical counterparts have 
probabilities of > 95% that they are associated with the X-ray 
source. We rejected objects that have probabilities of less than 
65%. All objects in the secure sample have probabilities >83%. 

After the X-ray counterpart verification our sample consists 
of 22 sources including 13 secure and 9 tentative objects. 
The separation between AGN and QSOs was made after the 
determination of the de-absorbed intrinsic X-ray luminosity (see 
Sect.0. 

Comments on the optical spectra of tentative objects: 
Marano 32A - all narrow emission lines have underlying broad 
(mainly blue-shifted) components. 

Marano 47A - the signal-to-noise ratio of the optical spectrum 
does not allow the exclusion of the presence of a weak, broad 
H/3 emission line. 

Marano 50A - H/3 and [Om] not covered by the spectrum, 
however, most likely a type II object since no Mg n but [O n] 
and [Ne v] emission lines. 

Marano 51 A - H/3 and [Om] covered by the spectrum but H/3 
coincides with the edge of the atmospheric A-Band. However, a 
strong broad H/3 emission line can be ruled out and [O n], [Ne v] 
emission lines are visible. 

Marano 66A - SNR of the optical spectra does not exclude the 
presence of a weak, broad H/3 emission line ([On] and [Nev] 
emission). 

Marano 116A - type II object/X-ray bright optical normal 
galaxy (XBONG) - galaxy spectrum with prominent [On] 
emission line, no Mgn, very weak [Om] emission lines, H/3 
coincides with the edge of the A-Band. 

Marano 133A - some narrow emission lines have underlying 
broad, blue-shifted components. 

Marano 253A - type II object/XBONG - no Mgn and H/3 
emission, weak [O m] emission, the second [O m] line falls into 
the atmospheric A-band. 

sdslb-014 - high SNR spectrum, H/3 and [Om] not covered 
by the spectrum, however most likely a type II object since no 
Mg ii but [On] and [Ne v] emission lines. 
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Fig. 1. R-band magnitude histogram of all 22 type II QSO can- 
didates (black solid line) and of the 13 secure type II QSO can- 
didates (grey filled histogram). 
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Fig.2. Observed 0.5-10keV X-ray flux vs. /?-band magni- 
tude. Large symbols represent secure type II QSO candidates 
while small symbols illustrate tentative type II QSO candidates. 
Dashed lines indicate fx/fopr values of 0.1, 1, and 10. 



3. Properties of the type II QSO candidate sample 

In TableQ]we summarise the observed properties of the objects. 
We list the name of the spectroscopically identified counterpart, 
optical coordinates, distance between spectroscopically identi- 
fied counterpart and X-ray source, the WFI-/?-band magnitude 
(unless otherwise mentioned), redshift, optical flag, X-ray coun- 
terpart probability, count rate, 0.5-10keV flux, and the Galactic 
absorption along the line of sight. 

The objects cover a redshift range of 0.5 < z < 3.278. 
Since the lack of emission lines in a wide optical wavelength 
range hampers the identification of type II objects, we only 
have one type II QSO candidate in the redshift interval of 
Z = 1-2. Figure [TJ shows the 7?-band magnitude histogram 
of the selected objects. As a simplification we consider the r- 
band (SDSS) magnitude to be equal to the 7?-band magnitude, 
although we are aware that shifts of 0.5 magnitude between 
R and r (SDSS) may occur. The majority of our sources have 
21 < R < 24. The 0.5-10keV flux in the sample ranges from 
0.3 -33 x 10- I4 erg/s/cm 2 . 



To calculate fx I f opt values we derived the optical fluxes in 
a band centred at 7000 A with a width o f 1000 A using the equa- 
tion / opt = i(T 04 *- 5 - 759 (|Zombecklll990h . X-ray fluxes were cal- 
culated for the 0.5-10keV energy range. The R-bwd magnitude 
vs. X-ray flux plane (Figure |2]i clearly shows AGN activity in 
all selected sources, since they have X-ra y- to-optical flu x ratios 
of /x//opt > 0.25. In accordance with lMainieri et alj {2002) 
we consider objects with fx/fopj > 0.1 as AGN. Almost half 
of the selected type II QSO candidates ha ve noticeable high X- 
ray-to-optical flux ratios (fx/fovv > 10). ISzokolv et all d2004) 
mentioned that type II AGN/QSOs cluster at higher X-ray-to- 
optical flux ratios than type I AGN. The majority of the spec- 
troscopically classified type II A GN/QSOs in the Lock man Hole 
also show high / x //opt values (Mainieri et al. 2002). The ob- 
served 0.5-10keV X-ray luminosities (not corrected for intrinsic 
absorption) of our objects range from L XoBS ~ 10 43 ~ 45 erg/s. 



4 M. Krumpe et al.: X-ray absorption in distant type II QSOs 



Table 1. Observed properties of type II QSO sample 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


No 


RA 


DEC 


'ox 


R 


z 


Opt. 


XID 


CR 


Flux 


Mt GAL 














flag 


probab. 




10~ 14 


10 22 








[arcsec] 










Iks-'] 


[erg s _I cm~ 2 ] 


[cm" 2 ] 


Marano 9A 


3 15 09.9 


-55 13 13 


1.18 


22.63 


1.427 


secure 


0.99 


21.3±0.9 


3.957 


0.027 


Marano 20A 


3 16 21.5 


-55 17 59 


0.63 


23.5 £ 


2.207 


secure 


0.97* 


17.1±1.3 


3.014 


0.027 


Marano 32A 


3 15 47.0 


-55 17 55 


1.44 


22.96 


2.727 


tentative 


0.99 


12.5±0.9 


3.497 


0.027 


Marano 39A 


3 13 39.7 


-55 01 51 


0.95 


23.41 


0.862 


secure 


0.98 


19.8+1.5 


3.982 


0.027 


Marano 47A 


3 15 38.7 


-55 10 44 


1.34 


23.94 


0.900 


tentative 


0.95 


4.5±0.5 


1.051 


0.027 


Marano 50A 


3 14 09.9 


-55 17 46 


1.25 


23.66 


0.986 


tentative 


0.96 


7.6±0.7 


1.667 


0.027 


Marano 51A 


3 16 30.6 


-55 15 03 


1.83 


20.95 


0.58 


tentative 


0.98 


10.4±1.0 


1.919 


0.027 


Marano 63A 


3 15 16.9 


-55 06 02 


1.25 


23.19 


2.800 


secure 


0.97 


7.2±0.7 


1.147 


0.027 


Marano 66A 


3 15 00.7 


-55 07 18 


1.23 


21.41 


0.981 


tentative 


0.99 


4.0±0.6 


0.566 


0.027 


Marano 116A 


3 16 20.7 


-55 16 52 


1.30 


21.2 £ 


0.581 


tentative 


0.99* 


5.4±0.9 


2.043 


0.027 


Marano 133 A 


3 14 26.2 


-55 21 13 


1.02 


23.54 


2.321 


tentative 


0.92 


2.7±0.6 


0.336 


0.027 


Marano 171 A 


3 13 51.4 


-55 02 56 


2.18 


22.17 


0.800 


secure 


0.86 


5.8±0.8 


1.998 


0.027 


Marano 224B 


3 13 04.8 


-55 16 04 


1.70 


21.51 


0.690 


secure 


0.96 


8.9+1.2 


3.899 


0.027 


Marano 25 3 A 


3 14 38.0 


-55 06 50 


2.64 


21.84 


0.517 


tentative 


0.65 


2.6±0.5 


1.341 


0.027 


Marano 463A 


3 16 25.3 


-55 08 39 


0.55 


25.0 £ 


2.531 


secure 


0.98* 


3.1+0.8 


0.789 


0.027 


Marano 61 OA 


3 15 51.8 


-55 12 22 


1.25 


21.12 


0.699 


secure 


0.98 


3.6±0.6 


1.275 


0.027 


X21516.135 


2 26 26.7 


-04 36 56 


0.90 


21.74 r 


3.278 


secure 


0.91 


0.9±0.2 


1.574 


0.025 


X0085 1.154 


22 15 41.6 


-17 37 54 


0.64 


19.00 


2.976 


secure 


0.96 


1.7±0.2 


1.210 


0.023 


X01 135.126 


1 52 57.6 


-14 08 40 


0.91 


20.87^ 


0.543 


secure 


g-jWFC 


10.5±0.5 


7.167 


0.015 


X03246.092 


43 46.0 


-20 29 56 


2.31 


20.01^ 


0.500 


secure 


0.83^° 


3.2+0.8 


4.104 


0.015 


phl5200-001 


22 28 26.4 


-05 18 20 


0.63 


22.0 r 


0.711 


secure 


0.99 WFC 


50.0+1.2 


33.03 


0.053 


sdslb-014 


2 18 43.0 


-05 04 37 


0.80 


22.0'" 


0.962 


tentative 


om wfc 


6.1+0.4 


6.380 


0.026 



Comments for Table [T] 

E - FORS R pre-imaging magnitude estimate; cross-calibrated with WFI-R-band magnitude catalogue; magnitude error estimate ±0.5 mag. 

K - since no R-band magnitude catalogue entry existed, the X-ray-to-optical counterpart probability calculation (see Eq. [T} is based on the S"-band magnitude 
catalogue. 

r - SDSS r-band magnitude (AB magnitude). 
UK - UK-Schmidt red plate magnitude. 

WFC - since no i?-band magnitude catalogue entry existed, the X-ray-to-optical counterpart probability calculation (see Eq.fTJ is based on the WFC-data. 
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4. XMM-Newton observations and extraction of 
X-ray spectra 

All data were pro cessed with the S AS version 7.0 (Science 
Analysis Software, Gabr iel et al.l2004l) package and correspond- 
ing calibration files. The epchain and emchain tasks were used 
for generating linearised event lists from the raw PN and MOS 
data. For all sources the effects of photon pile-up was negligible. 

The XMM-Newton data reduction for the M arano field sur- 
vey is described in detail in Krum peet alj d2007i) . For the XWAS 
and AXIS sources we downloaded all available X-ray data 
from the XMM-Newton archival] up to and including July 2007. 
Periods of high background were excluded from the analysis of 
all relevant data sets in the standard way. 

Circular or box-shaped source and background regions were 
manually determined for all contributing observations. Sources 
at large off-axis angles in the contributing observations were not 
considered as follows. The largest offaxis angles were 720 arcsec 
for the PN and 820 arcsec for the MOS detectors, respectively. 
The lower PN area was a result of the enhanced background con- 
tamination near the edges of the PN detector. Additional X-ray 
sources in the background regions were masked out. The auxil- 
iary response file (arf) was computed for each source and obser- 
vation individually. For the Marano field sources, we used ap- 
propriately "canned" response matrices from the XMM-Newton 
calibration homepag^ (XMM-revolution 110 and pattern 0-12). 
The PN y-coordinate of an X-ray source was determined to link 
the relevant PN response matrix file (rmf for single and dou- 
ble events, version 6.8) to the X-ray spectrum of the source. 
For the XWAS and AXIS sources the response matrix files were 
computed on a case-by-case basis. Where multiple X-ray spec- 
tra of a given source were added, the mean rmf was computed 
as a w eighted mean. For details of the procedure see Page et al.l 
(2003). The two MOS spectra were always added to form a sin- 
gle MOS spectrum. 

The X-ray source of the secure type II object phl5200-001 is 
surrounded by diffuse X-ray emission. For the reduction of the 
X-ray spectrum we used a smaller extraction radius for the point 
source and the diffuse X-ray emission as the background region. 



5. X-ray spectral analysis 

The X -ray spectral analysis was performed with XSPEC dArnaudl 
1996) version 12.3.0. Although we have a few objects with sev- 
eral hundreds of net PN source counts in the 0.2-8 keV range, the 
distribution peaks at ~40 net PN source counts. An appropriate 
X-ray spectral analysis for the low count regime has to be found. 
Although we only refer to the net PN counts in our simulations 
and in Table g] the fit uses both the PN and MOS data. Hence, 
the total number of counts used by the fit is typically twice that 
given in Table [4] 

5.1. Defining the appropriate fit statistic and binning method 

A problem notorious to X-ray astro nomy is pr o per fi tting 
of spectra with relatively few counts. iTozzi et ail d2006l) ap- 
proached this problem for the X-ray sources in the Chandra 
Deep Field South by running simulations for two different fitting 
procedures: Cash-statistics (unbinned) and classic ^-statistics 
with a binning of 10 counts per bin (min 10). They concluded 



1 http://xmm.esac.esa.int/xsa 

2 http://xmm.esac.esa.int/external/xmm_sw_cal/calib 



Table 2. Set of input parameters for the simulated X-ray spectra. 



Parameter 


value 


0.2-8 keV net counts 


10, 40, 100, 130, 200 


Redshift 


1, 2, 3 


Column density [cm -2 ] 


0, 10 21 , 10 22 , 10 23 , 10 24 


Photon index T 


2.0 



Table 3. Fit statistics and binning methods for the simulated X- 
ray spectra 3 . 



Fit statistic 


binning methods 


Cash-statistic 


unbinned, min 1, min 2, min 3, 




min 5, min 10, min 15, channel 60, 




channel 120, channel 180, channel 240 


^-statistic 


min 10, min 15 


Photon index T 


free fit, 2.0 



that the unbinned Cash-statistic fits recovered the input values 
better for X-ray spectra with less than 50 counts. 

We carried out a much more extensive study of the fit results 
as a function of fit statistic and binning method. This investiga- 
tion can be used to study the error distribution of the intrinsic 
column density (Nh) and to determine how many PN net counts 
are required to perform free fits of Nh and photon index F with 
acceptable errors in both parameters. 

We assumed emission from the AGN to be described as 
a power law with photon index F, modulated by Galactic 
foreground absorption and further i ntrinsic cold absorpti on at 
the redshift of the AGN. Following iMainieri etaTI d2002l) and 
Mate os et all (_2005), we simulated X-ray spectra with an input 
value of T = 2. All simulated X-ray spectra included Galactic 
absorption of Nh gal = 2 x 10 20 crrT 2 and considered all possible 
parameter combinations of Table [2] 

The X-ray spectra were normalised to reach the desired 
0.2-8 keV PN net counts with a deviation of up to 5%. We added 
Poisson noise to the X-ray spectra. We used a typical representa- 
tive background file (Marano 9A). As for the real data we added 
MOS 1 and MOS2 spectra to form a single MOS spectrum. The 
same source was simulated 1000 times for each set of specific 
parameters. 

The simulated X-ray spectra were grouped by using different 
binning methods (Table|3j. Cash and^ 2 -statistic were applied to 
recover the Nh input value. The redshift and the Galactic absorp- 
tion were set to the input value. The nominal initial guess for the 
intrinsic absorption was Nh = cm -2 , but other initial guesses 
up to Nh = 10 24 cirr 2 were tested as well. 

We studied the recovered Nh distribution with different fit 
statistics and binning methods by analysing the number of mis- 
matches, the peak, the shape and the significance of Nh detec- 
tion. The best fit statistic and binning method was selected by 
these criteria in the given order. In summary, the most important 
results of our simulations are as follows. 

- With only very few exceptions, the Cash-statistic with a bin- 
ning of one count per bin (min 1) recovered the input values 
best. 



3 The data were grouped from a minimum PHA channel (correspond- 
ing to 0.2 keV) to the maximum PHA channel (corresponding to 8 keV), 
with at least n counts per bin (abbreviated as min n) or with a fixed num- 
ber n of channels (channel n). 
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Even at a level of 40 net PN counts, the Cash-statistic 
with a binning of minimum one count per bin correctly 
recovers more than 90% of the objects with an intrin- 
sic absorption of A^h ~ 10 22 - 10 24 cm~ 2 up to z — 3.5. 
Figure [3] shows the retrieved A^h distribution for an in- 
put value of A^h = 10 22 crrT 2 at z = 2. The variance 
of the recovered A^h distribution increases with redshift. 
However, individual peaks clearly separate in the stud- 
ied redshift range, when the intrinsic A^h values differ 
by a factor of approximately 10 in hydrogen column 
density. The recovered A^h distribution for an intrinsi- 
cally unabsorbed source at z = 1 is shown in Fig. [4] 
Between 50-60% of all objects are correctly recovered 
independent of redshift. Although not all fits converge at 
A^h = cirT 2 , none of the resulting A^h values is signifi- 
cantly different from zero. The peak of the misclassified 



400 



300 



£ 200 



100 




absorbed objects shifts from Nu 
N H ~7x 10 21 cnr 2 forz = 3. 



10 cm forz = 1 to 



- For 10 net PN counts the fit retrieves a very broad A^h 
distribution for objects with an intrinsic absorption of 
Nh ~ 10 22 - 10 24 cirT 2 . Independent of redshift the fit 
recovers roughly 30% of the intrinsically absorbed X-ray 
sources as unabsorbed sources (Nh = cm -2 ). However, 
the fit does not significantly overpredict the A^h values. 
Although no absorption is recovered with a significance 
above 2<x due to the large errors, the fitted A^h values can 
be used as an estimate for the intrinsic A^h values. An un- 
absorbed X-ray source is correctly recovered by the fit in 
50-60% of all cases. In 90% of the cases the fitted N H val- 
ues for an intrinsically unabsorbed X-ray source do not 
exceed A^h ~ 10 22 cirT 2 forz = 1 andA^ H ~ 3xl0 22 cm -2 
for z = 2. 

- The unbinned Cash-statistic as used in iTozzi et al I (120061) 
shows the narrowest distribution in A^h for most of the pa- 
rameter combinations, but the peak of the recovered distribu- 
tion was found to be strongly dependent on the initial guess 

\12 



Fig. 3. Recovered A^h distribution for an X-ray source with 40 
net PN counts at z = 2. The input A^h = 10 22 crrT 2 is indicated 
by the dashed line. Fit method used: fixed F = 2, Cash-statistic, 
binning: min 1; 1000 simulations. 
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of JV H and weakly favours N H = 10 22 cm 2 independent of Fig 4 ReC overed N H distribution of an unabsorbed X-ray source 



input parameter for net PN counts less than 100 

- For a binning of min 10 or min 15 we find no difference 
between the Cash and the ^-statistic even at a level of 40 
net PN counts. Only the significance of the absorption is 
marginally higher for the Cash-statistic. There is also no dif- 
ference between a binning of min 10 and min 15. 

- Free fits in A^h and T are acceptable above 100 PN net counts. 
The recovered A^h -values show a narrow distribution at the 
input A^h value (FWHM~0.4dex). For all combinations of 
fitting and binning methods we retrieve the peak of T 2. 
For a binning of minimum one count per bin (min 1) the re- 
covered F distribution clearly peaks at T=2 with FWHM~0.8 
(FWHMnocounts ~0.4, FWHM 200 cou„ts -0.2). 



Simulations as shown in Fig. [3] and Fig. [4] are used to de- 
termine the A^h range that belongs to 68% and 90% of all 
simulations. In Fig. [5] we show these distributions for simula- 
tions at different redshifts (z=0,l,2,3,3.5) and column densities 
(Nh = 0, 10 22 , 10 23 cirT 2 ). The simulations contain 40 net PN 
counts and so are representative of many of our sources. The 
Cash-statistic with a binning of at least one count per bin is an 
appropriate method to determine the column density over a wide 
range of redshifts. It does not significantly overestimate the in- 
put A^h values. The fit recovers more than 90% of the unabsorbed 
sources up to z = 3 with A^h values below A^h = 10 22 cm" 2 . 



with 40 net PN counts at z - 1 ■ The input A^h = cm is in- 
dicated by the dashed line. A^h = cm 2 is represented by the 
lowest bin. Fit method used: fixed F = 2, Cash-statistic, binning: 
min 1; 1000 simulations. 



5.2. Modelling the X-ray spectra 

Following our analysis in Sect. I5.1l we grouped the extracted X- 
ray spectra in bins of at least one count per bin (0.2-8 keV) and 
used the Cash-statistic to determine A^h- The initial guess was 
set to A^h = 0. We corrected for Galactic absorption and per- 
formed free fits in A^h and F for spectra with more than 100 net 
PN counts. For X-ray spectra with less than 100 net PN counts 
we fixed F = 2 and only fit A^h- In these cases we give two errors 
for the value of A^h in Table [4] The first is the lcr error of A^h 
based on the fit with fixed F = 2, while the second takes into 
account the systematic shift for different photon indices (devia- 
tion in A^h for a fixed F — 1.7 and F = 2.3). The observed (not 
corrected for intrinsic absorption) and the intrinsic 0.5-10keV 
X-ray luminosities were computed by either a free fit in A^h and 
T for the bright sources (net PN counts >100) or a fixed F — 2 fit 
for the faint sources. 

Due to contamination by soft detector background we only 
used the 0.3-8 keV band for the X-ray spectral analysis of X-ray 
source X03246D92 (24 net PN counts in 0.3-8 keV band). The 
soft energies in the X-ray spectra of Marano 51A and Marano 
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Fig. 5. Diagnostic Nu vs. redshift plot. The plotted contours con- 
tain 68% (grey) and 90% (black) of the simulated X-ray sources 
with 40 net PN counts in the 0.2-8 keV band (NOTE: The fits 
used more than twice as many counts since both the PN and 
MOS data are used). The solid lines are the limits for intrinsi- 
cally unabsorbed sources. The limits for an intrinsic absorption 
of 7Y H = 10 22 
represent Nu 



i- 2 
10 23 



cm " are shown by dashed lines while dotted lines 



cm 



-2 



61 OA could not be well fitted by an single absorbed power law. 
In these cases a fit of an absorbed power law plus a soft excess 
component reproduce the X-ray data better. Therefore, we used 
a two component fit, an absorbed power law and an unabsorbed 
power law. No evidence for soft excess was found in the X-ray 
spectra of the brightest objects. The best fit models of the X-ray 
data are shown in Appendix lAl 

Table [4] shows the spectral parameters. In the columns we 
list the name of the spectroscopically identified counterpart (1), 
the XMMU source name (2) of the XMM-Newton X-ray source, 
the observed (3) and intrinsic (4) X-ray luminosity, the PN count 
number (5) in the 0.2-8 keV band, the absorbing hydrogen col- 
umn density (6), the photon index (7), and the quality for the 
used fit (8) when the X-ray data are modelled by a power law 
plus intrinsic absorption model, and finally, the photon index (9) 
and the quality for a purely reflection-dominated model (10). If 
no error is given for the photon index, this parameter was fixed 
for the fitting. 
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Table 4. Computed properties of type II QSO sample 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


No 


XMMU J 


l0 g( L X Bs) 


log(L XlNT ) 


-^Counts 


N s 


r 


C-stat/d.o.f 


r 


C-stat/d.o.f 






0.5-10keV 


0.5-10 keV 


PN 


10 

[cm -] 


intr. 
absor. 


intr. 
absor. 


pexrav 


pexrav 


Marano 9A 


031510.1-551313 


44.50 


44.55 


270 


0.3±0.1 


1.7+0.1 


426/503 


2.2±0.3 


473/503 


Marano 20A 


031621.6-551759 


44.82 


44.99 


64 


1.7±0.4±0.7 


2.0 


183/227 


2.0 


216/227 


Marano 32A 


031547.2-551755 


44.88 


45.22 


129 


9.8+2.4 


1.8±0.2 


285/291 


2.3±0.6 


295/291 


Marano 39A 


031339.7-550151 


44.04 


44.22 


36 


0.8±0.2±0.4 


2.0 


111/117 


2.0 


136/117 


Marano 47A 


031538.8-551043 


43.49 


43.72 


38 


1.5+0.5±0.5 


2.0 


149/157 


2.0 


180/157 


Marano 50A 


031410.1-551746 


43.76 


44.05 


50 


2.5±0.6±0.8 


2.0 


182/190 


2.0 


195/190 


Marano 51A 


031630.6-551501 


43.37 


43.61 


47 


2.5+0.9±0.9* 


2.0 


131/181 


2.0 


132/183 


Marano 63A 


031517.1-550602 


44.57 


44.81 


29 


4.1+1.4+1.4 


2.0 


125/140 


2.0 


127/140 


Marano 66A 


031500.8-550718 


43.36 


43.45 


39 


0.3±0.2±0.2 


2.0 


127/130 


2.0 


129/130 


Marano 116A 


031620.9-551651 


43.31 


43.72 


28 


5.6±1.2±1.0 


2.0 


129/126 


2.0 


136/126 


Marano 133 A 


031426.4-552113 


44.04 


44.04 


14 


0.0±0.8±0.0 


2.0 


87/70 


2.0 


94/70 


Marano 171 A 


031351.2-550257 


43.61 


44.00 


14 


5.4+1.5+1.0 


2.0 


46/51 


2.0 


63/51 


Marano 224B 


031304.9-551606 


43.71 


44.25 


8 


12.4+3.5+1.7 


2.0 


48/53 


2.0 


70/53 


Marano zjjA 


\Jj 14 ia.I- J 30O46 


42. y 1 


A1 CC 


1 1 
1 1 


14.j±j.U±j.z 


IX) 


04/00 


Z.U 


/ J/OJ 


Marano 463A 


031625.3-550839 


44.37 


44.53 


7 


1.6±1.0±0.9 


2.0 


76/45 


2.0 


81/45 


Marano 61 OA 


031552.0-551222 


43.28 


43.75 


17 


10.1+3.1+2.3* 


2.0 


110/119 


2.0 


113/120 


X21516.135 


022626.7-043654 


44.73 


45.15 


29 


14.8±5.7+4.8 


2.0 


61/79 


2.0 


66/79 


X0085 1.154 


221541.6-173753 


44.48 


44.70 


261 


5.6±2.0 


1.6±0.1 


572/621 


1.8±0.6 


583/621 


X01 135.126 


015257.5-140839 


43.85 


43.85 


392 


0.0+0.1 


1.7+0.1 


355/513 


2.1+0.2 


354/513 


X03246.092 


004345.8-202955 


43.48 


43.84 


35 


3.8+1.2+1.3 


2.0 


74/67 


2.0 


67/67 


phl5200-001 


222826.4-051820 


44.66 


44.92 


772 


4.5±0.4 


1.6±0.1 


834/970 


1.2±0.2 


1200/970 


sdslb-014 


021842.9-050437 


44.35 


44.45 


491 


0.4±0.1 


1.8+0.1 


546/593 


2.0±0.2 


650/593 



Comments for Table [4] 

* - a two component fit was used, an absorbed power law (Nu is given here) and an unabsorbed power law to model the soft excess 
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Fig. 6. Intrinsic A^h distribution of all 22 type II objects (black) 
and of the 14 type II QSOs (grey filled histogram). 



6. Discussion 

In Fig. [6] we show the fitted intrinsic jVh distribution of the 
sources. A column density peak at A^h = 4 x 10 22 cm" 2 is found. 
We find moderate absorption in the majority of our objects. The 
significance of the absorption exceeds 2cr in most of the cases 
(see the confidence contours in Appendix [A}. Two type II ob- 
jects are consistent with being unabsorbed X-ray sources, one of 
which is from the secure type II sample. 

The determination of the intrinsic A^h also allows us to 
compute the de-absorbed intrinsic X-ray luminosity Lx, NT °f 
our sources. As mentioned in Sect. [1] we define a type II ob- 
ject by the detection of narrow emission lines in the optical 
spectrum. Figure [7j shows the intrinsic column density vs. de- 
absorbed intrinsic X-ray luminosity plane. The dividing line of 
Lx iNT - 10 44 erg/s is used to distinguish between type II AGN 
and QSOs. Ten type II objects from the secure sample and four 
tentative objects are identified as QSOs. We detected one un- 
absorbed type II QSO but the object belongs to the tentative 
sample. No obvious trend of absorption in type II QSOs with 
intrinsic X-ray luminosity is found. 

Most of the type II QSOs fall into the same region 
of the Nh-L Xint diagr am where previous studies have also 
found type II QSOs dMainieri et al J 12 0021: ISzokolv et al.l l2004t 



lLa Franca et aLll2005t iPtak et al 



Mainieri et al 



2006). The additional criterion 



2002) makes only a small dif- 



(A^h > 10 22 cm" 2 
ference to our sample selection (two more secure objects). 

The A^H-redshift plane for type II QSOs and type II AGN is 
shown in Fig. [8] We found no obvious differences between the 
A^h distribution of AGN and QSOs. Although a tentative anticor- 
relation of Nu vs. redshift below z < 1 and a direct correlation 
above z = 1 is indicated, when all data is taken together there is 
no significant trend in A^h with redshift. Considering the typical 
flux and column densities found for our sources we expected all 
z > 1 type II objects to be classified as QSOs. At lower redshifts 
-40% of the type II QSO candidates are actually type II QSOs. 

Statistical fluctuations in the X-ray spectrum can lead to high 
values of spuriousl y measured A^h values at high redshifts (e.g. 
Akvl as et alj|200 6). However, Fig. [5] show that objects with in- 
trinsic A^h of several 10 22 - 10 23 cirT 2 are not significantly influ- 
enced by systematic trends in redshift. The fit method is also able 
to pick up much higher absorptions than found in our sources. 
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Fig. 7. Intrinsic A^h vs. de-absorbed intrinsic 0.5-10keV X-ray 
luminosity. Large symbols represent optical secure type II QSO 
candidates while small symbols illustrate the tentative sample. 
Open symbols indicate X-ray sources that have less than 40 net 
PN counts in the 0.2-8 keV band; filled symbols >40 net PN 
counts. The vertical solid line at log(Lx INX /(erg/s)) = 44 marks 
the dividing line between AGN and high luminosity QSOs. 
Objects with log(L XlNT /(erg/s)) > 44 and Nn > 10 22 cm~ 2 (up- 



per right corner) fall in the 
Mai nieri et al] d2002l) . 
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Fig. 8. Intrinsic A^h vs. redshift diagram for type II QSOs (top 
panel) and type II AGN (lower panel). Labels as in Fig. [7] The 
dashed and dotted lines represent the A^h values which corre- 
spond to a 30% and 50% decrease in the 0.5 keV flux at the given 
redshift, respectively. 



The scatter in the mentioned A^h range at z > 2.5 is consistent 
with an intrinsic absorption of A^h = 10 23 cm -2 for all sources. 

As a further test of the impact of statistical fluctuations, we 
plot in Fig. [8] the amount of A^h that is needed to reduce the 
0.5 keV flux by 30% and 50%. The 30%-line agrees well with 
the 90% contours of unabsorbed sources in Fig. [5] Consequently, 
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even for the sources at z > 2.5 the fitted A^h values are unlikely 
to be caused by statistical fluctuations. 

The previous conclusions depend on not having misinter- 
preted Compton-thick absorbed objects (A^h > 1.5 x 10 24 cm -2 ) 
as moderately absorbed objects. Increasing column densities 
cause a hardening of the X-ray spectrum. This trend stops when 
the material becomes optically thick to electron scattering. The 
X-ray spectrum of a Compton-thick absorbed source is com- 
pletely dominated by the reflection component. Their X-ray 
spectra show soft X-ray radi ation but with a much lower pho- 

ton in dex (T 4 - 1 3. iMaiolino et al.|[l998t iBassani et all 

Il999t iRisaliti et al.ll2000t Ewasawa et al.ll2001b IComasfrill2004h . 
We simulated an X-ray spectrum with F = 1.0 as observed 
in the 3-12 keV band for NGC1068. When we fit the spec- 
trum with a fixed T = 2.0, we recover values for absorption 
of N H = 0.7 x 10 22 cirT 2 , N H = 1.9 x 10 22 cirT 2 and Nn = 
3.8 x 10 22 cm~ 2 for redshifts of z = 1, 2, 3 respectively. Hence, 
moderately absorbed sources are found if intrinsically Compton- 
thick absorbed sources are studied. However, we are able to ex- 
clude a general misinterpretation of the spectra based on the fol- 
lowing arguments. 

First we tested if the X-ray spectra are better modelled with 
a purely reflection-dominated model (pexrav). We used solar 
abundances and a cutoff energy of 100 keV. The ratio between 
the reflected and direct component (fit parameter: scaling fac- 
tor) was set to 100 in order to obtain pure reflection. The cosine 
of the inclination was left as a free parameter. Again, we fitted 
the photon index as a free parameter only in the case the X-ray 
source had more than 100 PN net counts. 

Only the data of X-ray source X01 135.126 and X03246.092 
are slightly better represented by the reflection model in compar- 
ison to a single power law with intrinsic absorption. Further four 
X-ray sources (Marano 51 A, 63 A, 66 A, 61 OA) can be equally 
well modelled with a reflection model or a power law with in- 
trinsic absorption (see TableHJ). However, most of these objects 
have PN net counts of <50. Only X-ray source X01 135_126 has 
A^pN-Counts > 100. Low count spectra have the obvious prob- 
lem that their X-ray data can be well fitted by different models. 
Nevertheless, the majority of our objects is only adequately fitted 
by a model that includes intrinsic absorption. 

Secondly the reflected component of Compton-thick ab- 
sorbed sources is ~50-15 times weaker than the d e-absorbed in- 
trinsic X-ray luminosity ( Bra ndt & Has inger 20 05b. Hence the y 
are usually found at low values of fy/f[nm] dBassani et al.ll 999). 
For the majority of our objects the [O in] line is redshifted out of 
our spectral range. As a second best estimator we adopt fx/ f opt 
rather than /x//[oni]- Figure [2] shows that almost all objects 
have high /x//opt which is inconsistent with Compton-thick ab- 
sorbed objects. Only X0085 1_154 has a rather low /x//opt ratio, 
but the best X-ray spectral fit is a moderately absorbed power 
law with F = 1.6. Under the assumption of Compton-thick ab- 
sorption, all studied objects would be classified as QSOs and 
our sample would contain the most X-ray luminous QSOs ever 
studied (Lx INT ~ 10 47 erg/s). 



6. 1 . Stacking of X-ray spectra 

Another indication for Compton-thick absorption is the detec- 
tion of the iron line which is normally outshone by the direct 
component. Due to the reflection of the X-ray radiation from 
cold material in the torus, an iron Kor fluore scence line with an 
equivalent wi dth EW > 1 keV is expected dTurner et al J 1 1 9971 
Risaliti 2002). 
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Fig. 9. Stacked X-ray spectrum of all 14 type II QSOs (upper 
panel) and 8 type II AGN (lower panel). The QSOs stack con- 
tains ~1600 counts in the shown 4-8.2 keV energy range, the 
AGN stack 300 counts. The squares represent the ratio between 
the flux of the stacked spectrum and the simulated continuum 
spectrum. The solid line shows the lcr detection error while the 
dotted line represents the 2<x detection error. The solid grey line 
at a ratio of one illustrates the continuum flux. 



The low count numbers in our X-ray spectra do not allow the 
study of individual spectral features such as an iron line. Based 
on the assumption that the intrinsic X-ray properties in the sam- 
ple are similar, a stack of the individual spectra can reveal the 
presence of an iron line. Since the number of individual X-ray 
spectra increases the SNR in the stacked spectrum, we decided 
to include all 14 type II QSOs regardless of the optical flag. The 
assumption is justified because the secure and tentative sample 
do not show differences in fx/fopj (Fig.© and absorption val- 
ues (Fig^k 

ICorral et"ai] d2008) explain the stacking process in detail. 
Here we outline the essential components of the process. Using 
the best-fit parameters, unfolded spectra (eufspec in XSPEC) 
from the ungrouped observed X-ray spectra (MOS and PN data) 
are extracted. The unfolded spectra are corrected for Galactic 
foreground absorption, shifted to the rest-frame and normalised 
to the same rest-frame flux in the 2-5 keV band. The flux is re- 
binned into a common energy grid for all spectra. The stacked 
spectrum is binned to a final energy grid that contains at least 
100 counts per energy bin. The final errors of the stacked spec- 
trum are based on Gaussian propagation of the errors in the indi- 
vidual spectra. To distinguish between real spectral features and 
artifacts from the averaging process the underlying continuum 
has to be modelled. Each source was simulated 100 times us- 
ing the same model (absorption plus powerlaw) with the same 
2-8 ke V flux as observed in the real spectra. The individual sim- 
ulated X-ray continuum spectra from all different X-ray sources 
are stacked exactly as the observed X-ray spectra which are used 
to determine the lcr and 2<x contours (68% and 95% of all sim- 
ulated continuum spectra). 

The 0.5-10keV final stacked spectrum of the QSOs contains 
~4300 counts. Together with the lcr and 2cr contours the ratio 
of the observed spectrum to the average simulated continuum is 
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shown in Fig.|9j upper panel. The stacked type II QSO spectrum 
shows no residuals around 6-7 keV which would correspond to 
an iron Ka emission line. The SNR in the stacked spectrum is 
sufficient to detect an iron line with EW > 1 keV. The stacked 
spectrum does not provide any evidence for Compton-thick ab- 
sorbed sources based on a line emission and supports the as- 
sumption that our type II QSOs are actually moderately absorbed 
sources. 

Interestingly, the stacked spectrum of our 8 type II AGN does 
show a very prominent emission line (Fig. [9j lower panel), al- 
though in none of the individual AGN X-ray spectra a signif- 
icant excess at 6-7 keV is recognised (see the X-ray spectra in 
the Appendix [At. Since only a few objects are included in the 
AGN stack we used energy bins containing 30 counts instead of 
100 counts. We fitted the positive residuals in the stacked AGN 
spectrum with a Gaussian line profile. A line is detected with 
a significance above 2<x at a line energy of E - 6.6^'^ keV 
with a cr — 230^q° eV. The equivalent width of the fit is EW = 
keV. 

In principle there are two possible explanations for the de- 
tection of a broad iron fluorescence line. 

- a reflection line from the accretion disk, which can be signif- 
icantly broadened due to relativistic effects 

- a blend of narrow lines that is observed as a single unre- 
solved broad line due to the low signal-to-noise ratio of our 
observations 

The typical equivalent widths of the disk iron lines are usually 
significantly lower than the EW of the iron line found in our 
AGN stack. The energy and width of the line is more consis- 
tent with a blend of narrow lines. Depending on the properties 
of the reflecting material different ionisation stages of iron are 
excited as observed in the t he Compton-thick absorbed object 
NGC 1068 dOgle et al.ll2003l) . If the signal-to-noise ratio is low- 
ered the three resolved lines will be observed as an unresolved 
broad (EW > 1 keV) line with a peak above 6.4 keV. In this case 
a significant number of objects that we classified as AGN could 
be Compton-thick absorbed. Interestingly, except for one object, 
all X-ray sources, that can be better or equally well fitted by a 
purely reflection-dominated model in comparison to the power 
law plus intrinsic absorption model, are AGN. 

7. Conclusions 

We selected 5 1 spectroscopically classified type II objects from 
the XMM -Newton Marano field survey, XWAS, and AXIS. We 
re-investigate the optical spectra for narrow and common AGN 
high excitation lines, as well as verifying the X-ray counterpart 
determination. 22 sources meet our selection criteria (13 secure 
identification, 9 tentative). All sources have z > 0.5 and ob- 
served 0.5-10keV X-ray luminosities (not corrected for intrin- 
sic absorption) of Lx 0B s ~ 10 43 45 erg/s. The sample is not flux 
limited. 

The selection of both optical, narrow high excitation emis- 
sion lines and intrinsic X-ray luminosities ^Xint 

> 10 44 erg/s 

yielded 14 type II QSOs. 

Since the distribution of net PN counts peaks at ~40 the 
X-ray spectral analysis has to account for this very low num- 
bers of counts. We extensively simulated and studied such low 
count number X-ray spectra. A binning with at least one count 
per bin combined with the Cash-statistic recovered the input val- 
ues best. We proved that the method is able to find absorption up 
toN n = 10 24 cnr 2 . 



However, we discover only moderately absorbed type II 
QSOs. One QSO is consistent with being unabsorbed but it be- 
longs to the tentative sample. Compton-thick absorbed sources 
may be detected as moderately absorbed but we can exclude this 
scenario for the majority of our sources. The X-ray data are not 
well fitted by reflection models. The fx/ f opt values of the ob- 
jects and the non-detection of a broad iron line in the stacked 
type II QSOs spectrum give evidence that we have not misclas- 
sified Compton-thick absorbed sources. In contrary to the QSO 
stack, the stack of 8 type II AGN revealed a very prominent iron 
line with an EW ~ 2 keV. However, the shape of the single X-ray 
spectra and the modelling of the X-ray data showed that the ma- 
jority of our AGN are most probably moderately absorbed X-ray 
sources. 

The column density distribution found by us agrees 
well with those in deep Chandra and XMM-Newton surveys 
dMainieri et ail 120021; ISzokolv et all 120041: iMateos et all 120051: 
iPtak et all 120061) . These authors also reported a few cases of 
unabsorbed type II AGNs, as well as e vidence for some heav- 
ily absorbed sources (N H ~ 10 24 cm' 2 ). lLa Franca et all (120051) 
studied column density trends with X-ray luminosity in differ- 
ent redshift bins. They classified type II AGN as all sources 
that do not show any emission lines with FWHM > 2000 km/s. 
Independent of X-ray luminosity and redshift bins, they discov- 
ered type II AGN with an average absorption of A^h ~ 10 23 cm -2 , 
slightly above the column density peak of our survey (A^h = 
4 x 10 22 crrT 2 ). The present study of strictly classified type II 
AGN/QSOs, based on the optical spectra, verified that there is 
no obvious trend of absorbing column density with redshift or 
X-ray luminosity. 

O ur results apparently con tradict studies of th e local uni- 
verse, lassiniitai] (Q999|) and lRisaliti et all d2000l) found 75% 
of their sources with high absorption (A^h > 10 23 cm" 2 ) and 25- 
45% with A^h > 10 24 cirT 2 . Either the column density proper- 
ties change dramatically from the local to distant universe or the 
majority of the heavily absorbed distant sources are missed in 
our surveys. Indeed even the most luminous Compton-thick ab- 
sorbed sources in the local universe, e.g. NGC 6240, could not 
have been detected by our survey if they were at z > 0.4. 

In summary, the column densities in our survey show no 
trend in X-ray luminosity and no clear trend in redshift. If we 
compare our results with samples of AGN in the local universe, 
our sample does not contain a significant fraction of heavily ab- 
sorbed sources (A^h > 10 23 cirT 2 ). Our survey of objects with 
z > 0.5 is limited to observed X-ray luminosities in excess 
of Lxobs > 10 43 erg/s. We can only expect to find QSOs in- 
trinsically more luminous than Lx, MT > 10 45 e rg/s, if they are 
Compton-thick absorbed. Norm an et al.l ([2002) claim that for 
the most distant type II QSO ever detected (z - 3.7), there is 
strong evidence for heavy or even Compton-thick absorption. 
However, our survey rules out large numbers of Compton-thick 
absorbed sources with X-ray luminosities of Z>x INT > 10 45 erg/s. 
Hence, potential Compton-thick absorbed objects at high red- 
shifts are likely to have similar X-ray luminosities to Compton- 
thick absorbed objects in the local universe. In order to find a 
supposed, rare population of very luminous, Compton-thick ab- 
sorbed QSOs_jJarge£_survey area is needed. The 2XMM cata- 
logue dWatson e t al. 2008) with a survey area of ~360deg 2 could 
provide a valuable source to reveal such a population. 
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Appendix A: Optical & X-ray spectra, confidence contours 

In this section we show the optical and X-ray spectra, as well as the confidence contours of X-ray spectral fits for the parameters 
Af H andr. 

Optical spectra: Optical atmospheric absorption corrected, wavelength and flux calibrated spectrum for the optical X-ray 
counterpart is shown. All spectra are in flux units of 1CT 18 erg ctrT 2 s _1 A -1 . The black solid line represents the spectrum, the green 
solid line the error spectrum (not available in all spectra). Red markers indicate possible emission lines, blue markers absorption 
lines. There are exceptions in a few spectra. Spectral features at 5580 A are spurious due to incomplete subtraction of a night sky line. 

X-ray spectra: 

The X-ray spectra show the PN (black) and combined MOS (red) data, as well as the best fit model (foreground Galactic 
absorption, power law with intrinsic absorption at the object's redshift). In the case of Marano 51A and Marano 610A the shown 
best fit model consists of foreground Galactic absorption, an unabsorbed power law, and an power law with intrinsic absorption at 
the object's redshift to account for the soft excess in the X-ray data. 

For objects that have less than 100 PN source counts in the 0.2-8 keV band we fixed the photon index (T = 2). Otherwise, free 
fits in jVh and F are shown. The fit parameters are given in Table [4] For illustration purposes the X-ray data have been rebinned to 
different signal-to-noise ratio, after grouping to a minimum of one count per bin (min 1). 

Confidence contours: Confidence contours of the absorbing column density in units of Nh = 10 22 crrT 2 vs. photon index Y of 
the X-ray spectral fits. Confidence contours are plotted for 68% (black), 90% (red) and 99% (green). The contours are based on 
free A^h and F X-ray spectral fits for all objects independent of PN net counts. 

Comments: Below every set of optical & X-ray spectra and X-ray contour a comment for the objects shown is given. This 
comment includes the classification number of the optical counterpart and the corresponding redshift. Furthermore, a short fit 
description is given. For more details on the properties of the sources and the X-ray spectral fits see Tables Q~|[4] 
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object Marano 9A, z = 1.427, free X-ray spectrum fit 
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object Marano 20A, z = 2.207, frozen F = 2.0 X-ray spectrum fit 
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object Marano 39A, z = 0.862, frozen Y = 2.0 X-ray spectrum fit 
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object Marano 47 A, z = 0.900, frozen Y = 2.0 X-ray spectrum fit 
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object Marano 50A, z = 0.986, frozen F = 2.0 X-ray spectrum fit 
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object Marano 51 A, z = 0.58, two power laws, frozen F = 2.0 X-ray spectrum fit 
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object Marano 66A, z = 0.981, frozen T = 2.0 X-ray spectrum fit 
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object Marano 133A, z = 2.321, frozen F = 2.0 X-ray spectrum fit 



M. Krumpe et al.: X-ray absorption in distant type II QSOs, Online Material p 8 




5000 6000 7000 

WAVELENGTH (A) 





Energy (keV) Parameter: nH (10 22 > 

object Marano 171 A, z = 0.800, frozen F = 2.0 X-ray spectrum fit 
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object Marano 224B, z = 0.690, frozen F = 2.0 X-ray spectrum fit 
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object Marano 463A, z = 2.531, frozen F = 2.0 X-ray spectrum fit 
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object X00851_154, z = 2.976, free X-ray spectrum fit 



2.5 
2 



J III I 



i ii ( 



5000 6000 
WAVELENGTH (A) 





Energy (keV) Parameter: nH (10 22 > 

object X01 135J26, z = 0.543, free X-ray spectrum fit 



M. Krumpe et al.: X-ray absorption in distant type II QSOs, Online Material p 12 



fJ||]P|ffl ,f f l Jtt^, 



lull 



kkiwk 



5000 6000 7000 

WAVELENGTH (A) 




Parameter: nH (10 22 i 



Energy (keV) 

object X03246.092, z = 0.500, frozen T = 2.0 X-ray spectrum fit (0.3-8 keV) 
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object ph!5200-001, z = 0.711, free X-ray spectrum fit 
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object sdslb-014, z = 0.962, free X-ray spectrum fit 
(optical spectrum taken from Akiyama et al. 2008) 



